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Abstract 

The coordinate- momentum double distribution function p(r, p)d 3 rd 3 p is 
calculated in the local density approximation for bosons with positive scat- 
tering length a in a trap. The calculation is valid to the first order of a. To 
clarify the meaning of the result, it is compared for a special case with the 
double distribution function p w d 3 rd 3 p of Wigner. 



Using the local density approximation (LDA) [1,2], which is a straightforward adap- 
tation of the Thomas- Fermi method, the density distribution p(r)d 3 r in coordinate space 
for BEC for a > in a trap has been obtained. We want to calculate in this paper the 
coordinate-momentum distribution p(r, p)d 3 rd 3 p in the same approximation. We follow 
the notation of Ref. [2] throughout. In particular, the fugacity z of the system is 

z = exp[fi/kT] (1) 

where fi is the chemical potential. We introduce a local fugacity ((r) defined as 

C = zex P [-pV(v)]. (2) 

1. The Gaseous Phase 

By the gaseous phase we include both the system before BEC sets in, and the 
system at high densities for the cells outside of r [2], i.e., outside of the region where 
BEC takes place. We consider such a cell of volume V in which the local fugacity is (. 
Using the method of Ref. [3] we write the grand partition function in the cell as 

Q = J2C N Tr[exp(-^H - (3H% (3) 

N 

The average occupation number <C n k 3> of the state with momentum hk can be com- 
puted from 

Q « n k »= J2( N Tr[(ala k )exp(-(3H - f3H')}. (4) 

N 

We shall drop all terms beyond the first order of H' . Since H commutes with a\ak, we 
find 

Q = Y, ( N Tr[eM-PH )(l - (3H')\ = Q + Qi (5) 

N 

and 

Q « n k »= £ ( N Tr[exp(-(3H )(ala k )(l - f3H')\ =A + Ai (6) 

N 

where 

Qo = E ( N Tr[eM-PH )} =i[[l- (e^f 1 (7) 

N 



2 



is the term in Q without the perturbation term. Qi has been evaluated in Ref. [3]. 



Qi/Qo = -P 



mV 



Y n « n /3 + Y \ n l ~ 2 



where the bar means average over the grand canonical ensemble Qq: 
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Ao/ Qo = n k 



AJQo = -(3 
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where the symbol ( ) means the same average as the bar. The coefficient in (8) and (11), 
—f3Anah 2 (mV)~ 1 , is equal to — 2a\ 2 V~ 1 . Now (11) can be rewritten as 



A 1 /Qo = -2a\ 2 V- 1 ((n k -n k ) 



Y n « n /3 + Y \ n l ~ Y \ n <* 

a z a z 



+ n k Q 1 /Q . (12) 



Notice that Qo is a product distribution function according to (7). Thus n a rip = n a np 
if a ^ (5. Using this and similar identities we find that in the sum over a and (5 in (12), 
the bracket ( ) vanishes unless k = a or k = (3. Thus 

Ai/Qo = -2aAV _1 ( Y 2n k n p {n k - n k ) + \{n\- n h n\ -n 2 k + n 2 ) ) + n k Q 1 /Q - (13) 



(14) 



Now V 1 (E/^k 71*71/8(71* - n k )) -> p(n| - n£) as V -> oo, yielding 



Ai/Qo = -Aa\ 2 p(n 2 k - n%) + n k Qi/Q - 

Adding this to (10) and dividing by 1 + Q±/Qo we obtain, to order a, 

< n k >= n fc - Aa\ 2 p(n 2 k - 7j|). 

The number of modes k in a cell of volume V is (8ir 3 ) 1 Vd 3 k. Thus the combined 
coordinate-momentum distribution p(r, p) is given by 



(15) 



h 3 p(r, p) =« n k »= - - 4 7 rAV(r) (1 f e Ce / !; £a)2 



(16) 



where e k = and ( is given by (2). In (16) we have evaluated nf in a straightforward 
way from the product partition function (7). 

Integrating (16) over d 3 p we should get the density p(r) times h 3 . This can be done 
without much difficulty yielding Eq. (3) of Ref. [2]. 

2. The Region with Condensate 

For high densities, BEC forms in some cells of the trap. In those cells p = po + p s > 
p , where [4], 

Po = A- 3 £3/ 2 (l) (17) 

and 

V(r) + Anap s (r)h 2 /m = V{r ). (18) 

Here p s denotes superfluid density, i.e., density of particles with p = 0. An important 
parameter £ 5 = p s /p, a function of the location of the cell, with value between and 1, 
describes incomplete occupation of the ground state, and was studied in detail in Ref. [5]. 
[Notice that £5 and £ are totally different quantities.] For cells without BEC, £5 = 0. 

It was shown in Ref. [5] that the system in a cell with BEC has an energy given 
by (5.16) with a phonon spectrum (for k 7^ 0) given by (5.18): 

Ti 2 

hou k = — (k 4 + 2klk 2 ) 1/2 , k\ = 8vra£ 5 p = 8nap s . (19) 

Notice that for the gaseous phase, £5 = and the phonon spectrum is quadratic for small 
k. 

The phonon creation operator b k and the particle creation operator a\ are related 
to each other through a Bogoliubov transformation [6]: 

a k = (b k -a k bl k )(l-al)' 1/2 (20) 

where 

ot k = kv 2 (k 2 + k 2 - ^k A + 2k 2 k 2 ). (21) 



4 



For a state with rrik phonons we can compute the average occupation number <C n& ^> of 
atoms in the state p = Kk using (20) above. The result is linear in m^. Now the average 
number of rrik is given by Eqs.(5.27) and (5.31). Thus 

p(r, P ) = h~ 3 \a\ + (1 + ai)(e?**> - 1)^](1 - a 2 ,)' 1 , (k + 0), (22) 

where u>k is given by (19), and ak is given by (21). 

For k ^> k = y/8nap s , the phonon energy (19) can be expanded in powers of a 
and (22) becomes 

p(r,p) = /*-V^-ir\ (*;»*<,), (23) 

where 

For other values of k > 0, Eq.(22) gives the distribution. It is a complicated function of 
k. For < k <C k , it reduces to 

p(r, p) h^m^P' 2 , (0 < fc < fcb). (24) 

Notice that this differs by a factor of 2 from the corresponding distribution when a = 0. 

3. Wigner Double Distribution 

What is the meaning of the double distribution p(r, p)? It, of course, should only 
be used [2] for d 3 r > (L 2 ) 3 , and for d 3 rd 3 p > h 3 . But does it have a clear meaning in 
quantum mechanics? We discuss this by examining Eq.(16) in the limit of a = 0, for 
the case of a spherically symmetrical harmonic trap V(r) = ^muj 2 r 2 . In such a case we 



can compute exactly the matrix element of (r'| ^ ~"h l r ) — X^i^'l^e -18 ^ |r), by using, 
e.g., the result of Ref.[9]. Using Wigner's idea [10], we put r' = R— | rj, and r = R+ | r\ 
and evaluate the above, and then make a Fourier transform to the variable P conjugate 
to rj. 



The resultant double distribution a £a Wigner becomes 

-)=exp - (tanA )( 



p.(R,P) = /^/(sec/^) 3 exp -^(tan/ i ^)(— + -mu 2 R 2 ) (25) 



where e = (3Tiw. In the limit that e — > 0, this is exactly Eq.(16) for a = 0, [noticing that 
the local fugacity ( is given by (2)] which is in agreement with the discussion in Ref. [2] 
for the single distribution function p(r). 

The work of CNY is supported in part by an NSF Grant PHY-9309888. That of 
LHY is performed under the auspices of US DOE. 
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